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Abstract 

In this paper a method previously employed by Salmhofer to establish an 
exact equivalence of the one-flavour strong coupling lattice Schwinger model 
with Wilson fermions to some 8-vertex model is applied to the case with 
two flavours. As this method is fairly general and can be applied to strong 
coupling QED and purely fermionic models with any (sufficiently small) num- 
ber of Wilson fermions in any dimension the purpose of the present study 
is mainly a methodical one in order to gain some further experience with it. 
In the paper the vertex model equivalent to the two-flavour strong coupling 
lattice Schwinger model with Wilson fermions is found. It turns out to be 
some modified 3-state 20-vertex model on the square lattice, which can also 
be understood as a regular 6-state vertex model. In analogy with the one- 
flavour case, this model can be viewed as some loop model. 
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The motivation for the present work derives from a paper by Salmhofer [T|] 
who demonstrated the exact equivalence of the one-flavour strong coupling lattice 
Schwinger model with Wilson fermions to a 8- vertex model which can also be under- 
stood as a certain self-avoiding loop model. While strongly coupled gauge theories 
with staggered fermions lead to pure monomer-dimer systems || the considera- 
tion of Wilson fermions results in more complicated models, as the investigation 
of Salmhofer Q demonstrates by example. The particular equivalence established 
by him has led to hitherto unfeasible computer studies @-|| as well as to the 
application of certain analytical methods H yielding previously unknown informa- 
tion about the critical behaviour of the one-flavour strong coupling lattice Schwinger 
model. However, the method employed by Salmhofer is not specific to this model. 
It can equally be applied, under certain technical complications, to strong coupling 
lattice QED (in fact, to any gauge theory, with some modifications) and purely 
fermionic models with any (sufficiently small) number of Wilson fermions in any 
dimension. It is therefore particularly interesting for models with an odd number of 
fermion species which give rise to serious difficulties for any numerical simulation so 
far. Furthermore, the vertex model equivalence (which has to be found case by case) 
of certain lattice models of quantum field theory with Wilson fermions possibly will 
allow to make use of cluster algorithms for their numerical study in order to reduce 
the computational effort required. Beyond this, the equivalence may also lead to 
new ideas for analytical investigations. 

In this paper we extend the analysis of Salmhofer fl|] to the two-flavour strong 
coupling lattice Schwinger model. Beyond the specific result of the paper - the vertex 
model equivalence for the model under consideration - the main purpose of it is to 
collect some practical experience in dealing with models which have a larger number 
of Grassmann variables per lattice site than the one-flavour Schwinger model with 
Wilson fermions, which has four. The main difficulty is an algebraic one. While 
the study of Salmhofer requires paper and pencil only, to go beyond becomes 
possible only by using an appropriate formula manipulation program. In this study 
we are using a program for Mathematica JF] which can efficiently handle Grassmann 
variables and which has been designed just for the present purpose. 

Before going into technical details, let us shortly describe the method to be used 
in words ||T|. The square lattice A the two-flavour Schwinger model is defined on is 
divided into an even and an odd sublattice. Use will be made of the fact that effec- 
tively only fermionic degrees of freedom are left in the theory under consideration (In 
the strong coupling Schwinger model the bosonic gauge field degrees of freedom are 
integrated out exactly first, leaving the fermionic sector for further consideration.). 
First, on the (say) even sublattice A e the Grassmann integration is carried out. This 
is a local procedure (the reason for it being, that at most nearest-neighbour cou- 
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plings are present in the action) which amounts to a completely algebraic problem 
due to the Grassmannian character of the integration. The remaining Grassmann 
variables live on the odd sublattice A G and certain graphical rules, i.e. vertices, can 
be assigned to their combinations in a fairly natural way. In the final step an anal- 
ysis can then be made what structures (lattice clusters) built of these vertices the 
remaining Grassmannian integration on the odd sublattice leads to. Since for Wil- 
son fermions the lattice is homogeneous the result can be described in terms of a 
vertex model. 

We are now prepared to pursue the program for the two-flavour Schwinger model 
with Wilson fermions (with Wilson parameter r = 1). For easier reference, our 
notation follows closely that applied in M. The partition function Z\ of the model 
on the square lattice A is given by 



Z A = I DUUD^D^ e~ s , (1) 

i=l 

where DipiDipi = FLeA Ila=i dip iy0l (x) dtpi, a (x) denotes the multiple Grassmann in- 
tegration on the lattice (i is the flavour index.). The action S is defined by 



S = S F + PS G , (2) 
2 (\ 

Sf = EE uEfe+^Ki+^WiW 

xeAi=l \ Z A* 

+ $i(x)(l - ^)Ul(x)iPi(x + ej) - M^iix^iix)^ , (3) 

and Ufj, = exp [— iA^], (3 = 1/g 2 . Sq is the standard Wilson action and the hopping 
parameters Ki are given by K{ = l/2Mj (For the sake of generality we allow for 
different masses for each flavour here.). In the strong (infinite) coupling limit /3 — 
the gauge field integration can be performed exactly and one finds (cf. eq. (5) in JL]) 



£ ff[D^D^ exp(5:^M^(x)^(x)] 

fc ; e{0, 1,2,3,4} i=l \xeAj=l ) 



1 



x n 



k, 



£ Mx)T^ m {x + e M ) ) ( £ Mx + e,)fi(,) . (4) 



\m=l 



\n=l 
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In the two-flavour Schwinger model, the occupation numbers k\ cannot exceed 4 
because on each lattice site only 8 Grassmann variables are present. In fact, explicit 
calculation shows that by virtue of the nilpotency of the Grassmann variables fc/ 
cannot exceed 2. As in [|IJ, we choose 71 = <7 3 , 72 = <J\ (<jj are Pauli matrices.) and 
the projection operators Tjf> = (l + ej /J )/2, e 6 { — 1, 1} then read T± = diag{l, 0}, 

t[ ^ = diag{0, 1}, T 2 = iQi)- ^ * s now advantageous in order to simplify further 
results to define for each flavour combinations of the original ip fields in which terms 
is diagonal (The index in front of the comma always is the flavour index, the 
one behind it the component index.) 

1 T 

Xi{x) = Uipi(x) = -j= {^i,i{x) +ip ii2 {x),^i >1 (x) - ipi,2{x)) , (5) 



Xi{x) = ipi{x)U = (i>i,i{x) +ipi,2(x),i>i,i(x) -il) ii2 (xf) (6) 

with 




u = -e ; ; = u- 1 ■ (7) 



In the expression for the partition function (|4]) the integration on the even sublat- 
tice A e is now performed first. As each lattice site x supports 8 Grassmann variables 
terms which respect 



E h + s * = 4 (8) 

can give a non- vanishing contribution only (s x is the power of X^ =1 M^i^ipi^x) in 
eq. (§)•)• Explicit calculation now reveals that only terms with s x = or s x = 4 
are non-vanishing. A given set of k\ on the lattice A uniquely determines a certain 
configuration contributing to the partition function. Below we give the results for 
all local integrals related to a given lattice site x G A e which allow non-vanishing 
contributions. The results have been arranged in a way most suitable for the fur- 
ther discussion. While in the calculation for the one-flavour model has been 
performed by means of paper and pencil, the two-flavour case requires use of the 
computer and we have relied on a purpose-written Mathematica program to do 
the necessary algebra. The graphical rules below have to be interpreted as follows: 
the black dot denotes any point x on the even sublattice A e at which the Grass- 
mann integration is performed, a dashed line means k\ = while thin and thick 
lines stand for k\ = 1 and k\ = 2 respectively. The first coordinate component of 
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x is understood to be the horizontal one while the second component runs ver 



There is one possible vertex with s x = 4: 



Vertex 1 



- + - 



2 2 



!=la=l 4 - \j=l / 



Ml Mj 



All further vertices have s x = of course. 



Vertex 2 



t 



Vertex 3 



Vertex 4 



Xi,i(x + e 2 ) X2,2( x - e s) XiA x ~ e 2) X2,l( a; + e 2) 
x xi,2(a: + e 2 ) X2,i(^ - e 2 ) - e 2 ) X2,2(a: + e 2 ) 



^1,1 ( x + e ^2,2(2; - ei) ^1,2(2; - ei) ^2,1^ + e i) 
x ^i j2 (x + ei) ^2,1 - e ^1,1 - ei) ^2,2(2 + ei) 



r 



Xi, 2 (a; - e 2 ) $2,1 fa + e i) ^i,ifa + e X2,2(^ - e 2 ) 
x Xi.ifa - e 2 ) ^2,2(3; + ei) ^i )2 fa + ei) X2,ifa - e 2 )/4 



Vertex 5 
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Xi,i( x + e s) ^2,2(2; - ei) i/)i, 2 (x - ei) X2,ifa + e 2) 
x Xi^fa + e 2 ) ^2,1 (x - ei) ^l.ifa - e t ) x 2j2 (x + e 2 )/4 



Xl,l( X + e 2) 1p2,l( X + e *Pl,l( X + e X2,l(^ + e 2) 

x Xi,2(a? + e 2 ) ^2,2(3; + ei) ^1,2(3? + ei) X2,2(x + e 2 )/4 



Xi,2(a? - e 2 ) ^2,2(3: - ei) ^1,2(2: - ei) X2,2(x - e 2 ) 
x Xi,i(^ - e 2 ) ^2,1 (x - ei) ^1,1 (x - ei) X2,i(x - e 2 )/4 



x 



x ip 



Xl,l(x + e 2 ) X2,2( X ~ e 2) 

Xi,2(x + e 2 ) X2,i(x - e 2 ) + xi,i( x - e 2) X2,2(x + e 2 ) 
''1,2(2! - ei) ^2,2(2: - ei) - e i) ^2,1 - ei)/4 



+ XiA x - e 2) X2,i(x + e 2 ) 



X 



X *0 



Xi,i(x + e 2 ) X2,2(x - e 2 ) + xiA x ~ e 2) X2,i(x + e 2 ) 
Xi,2(x + e 2 ) X2,i(x - e 2 ) + xi,i(^ - e 2 ) X2,2(x + e 2 ) 



•i,i(x + ei) ^2,i(x + ei) ^1,2(3? + ei) ^2,2^ + ei)/4 



^i,i(x + ei) ^2,2(3? -ei) + ^1,2(3? - ei) ^2,i(x + ei) 

x |Vi,2(3: + ei) ^2,1(3: - ei) + Vi.ifa - e ^2,2(3: + ei) 
x + e 2) X2,i(x + e 2 ) xi,2(x + e 2 ) X2,2(x + e 2 )/4 



i/)i,i(x + ei) ^2,2(3; - ei) + ^1,2(3; - ei) ^ 2) i(x + ei) 
x [^1,2(3: + ei) ^2,1(3: - ei) + Vi.iC^ - e ^2,2(3; + ei) 
x XiA x - e 2) X2,2(x - e 2 ) xi,i(a; - e 2 ) x 2 ,i(z - e 2 )/4 
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Vertex 12 



Vertex 13 



XiA x - e ?) i>2,i( x + e + ^1,1 + e i) X2,2(a; - e 2 ) 
x |xi,i(x - e 2 ) ^2,2^ + ei) + ipi, 2 (x + ei) X2,i(a? - e 2 ) 
x ^i, 2 (a:-ei) ^2 >2 (x-ei) ^i,i(x-ei) ^2,i(x-ei)/2 (20) 



Xi,2(a? - e 2 ) V2,i(^ + ei) + + ei) X 2)2 (x - e 2 ) 

- e 2 ) ^2,2^ + ei) + ^1,2^ + ei) x 2 ,i(x - e 2 ) 
x + e 2 ) X2,i(x + e 2 ) Xi,2^ + e 2 ) X2,2(x + e 2 )/2 (21) 



Vertex 14 



Vertex 15 



L Xi,i(z + e 2 ) ^ 2)2 (a; - e x ) + ^1,2^ - ei) x 2 ,i(x + e 2 ) 

x Xi,2(z + e 2 ) ^2,1 (x - ei) + Vi.ifa - e X2,2(a; + e 2 ) 
x ^1,1 (x + ei) ^2,1 (x + ei) ^1,2(2; + ei) ^2,2(3: + ei)/2 



(22) 



Vertex 16 



Xi,i(x + e 2) ^2,2(2; - ei) + Vi,2(x - ei) x 2 ,i( a; + e 2) 
x [xi,2(x + e 2 ) ^2,1 (a; - ei) + ^i,i{x - e x ) X2,2(x + e 2 ) 
x Xi,2(x - e 2 ) X2,2{x - e 2 ) xi,i(x - e 2 ) X2,i(x - e 2 )/2 



(23) 



^i,i(x + ei) X2,l( a; + e 2) + Xi.i^ + e 2) i>2,i( x + e i) 
x [^i,2(x + ei) X2,2(x + e 2 ) + Xi,2(x + e 2 ) ^2,2^ + e x ) 
x ^1,2 (x - ei) ^2,2(2; - ei) ^1,1 - e i) ^2,1 (x - ei)/2 



(24) 
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Vertex 17 



Vertex 18 



i>i,i( x + e X2,i( x + e 2 ) + + e 2 ) ^2,1(3: + ei) 

L ^1,2(3; + ei) X2,2(a: + e 2 ) + Xi^tx + e 2 ) ^ 2 , 2 (a; + e x ) 
x xi,2(a: - e 2 ) ^2,2(2; - e 2 ) xi,!^ ~ e 2) X2,i(z - e 2 )/2 



Vertex 19 



Vertex 20 



^1,2(3? - ei) X2,2(x - e 2 ) + Xi,2(a? - e 2 ) ^2,2(3? - ei) 

"01,1 - e X2,i( a; - e 2) + Xi.iC^ ~~ e 2) ^2,1 - ei) 
x xiA x + e 2) X2,i(x + e 2 ) xi,2(a: + e 2 ) ^2,2(3; + e 2 )/2 



$i, 2 (x - ei) x 2)2 (a; - e 2 ) + xi,2^ - e 2 ) ^2,2(3; - ei) 

x ^1,1 _ e X2,i( ;r - e 2) + Xi.iC^ ~~ e 2) ^2,1(2; - ei) 
x ^1,1 (x + ei) ^> 2 ,i(x + ei) ^>i, 2 (x + ei) ^2,2(3; + ei)/2 



x 



X 



X 



X 



X 



Xi,i( x + £2) i>2,2{x - ei 

Xi,2(^ - e 2 ) ^2,1 + e i 
^i >2 (a; - ex) x 2j2 (:r - e 2 

^i,i(x + ei) X2,i(x + e 2 

Xi,2(a! + e 2 ) ^2,1 (x - ei 

Xi,i(^ - e 2 ) ^2,2(3; + ei 
^1,1 (x - ei) X2,i(x - e 2 

^1,2(3? + ei) X2,2(x + e 2 



+ ^1,2(3: - ei) X2,i(^ + e 2 ) 

+ ^1,1(3? + ei) X2,2(x - e 2 ) 

+ Xi,2(x - e 2 ) ^2,2(2: - ex) 

+ X\,\{ x + e 2) ^2,1(3: + ei) 

+ ^1,1 (a; - ei) x 2 ,2(z + e 2 ) 

+ ^1,2(3; + ei) X2,i(x - e 2 ) 

+ - e 2 ) ^2,1 (a: - e\) 

+ Xi,2(x + e 2 ) ip 2 ,2( x + ei) 



Remarkably, as one recognizes from the above equations the remaining field combi- 
nations on the odd sublattice factorize for each vertex into exactly identical contri- 
butions from the ip, X and the i/j, x subsystems. To see this use the map 



entailing 



< * i>i,2(x) 

Xi,i{ x ) < — ► XiA x ) > 
XiA x ) < — ► — Xi,i(x) 



(29) 



(30) 



The above mentioned factorization property for each vertex comes as a surprise be- 
cause it could hardly have been guessed prior to the explicit calculation. The vertex 
factorization property immediately leads to the conclusion that any individual con- 
tribution (of any vertex cluster) to the partition function Z\ is positive. The minus 
sign in front of the r.h.s. of eqs. flnp-flTUp, (B^)-(B^) at first glance disturbing this 
property is properly taken care of by the map ([29]) , (|30|) . It should be mentioned 
that a similar vertex factorization property holds in the one-flavour Schwinger model 
jjH, there however being fulfilled rather trivially. One might wonder, whether a ver- 
tex factorization property also holds for the strong coupling Schwinger model with a 
larger number of fermion species than two, perhaps being a rather general structure. 

The vertex factorization property means that we can now restrict ourselves to 
the discussion of the tp, x subsystem. Each vertex cluster contributing to the parti- 
tion function Z\ adds a term to the partition sum equal to the square of the weight 
of the corresponding cluster in the ip, x subsystem. We consequently need to know 
weights in the ip, x subsystem up to minus signs only. In order to achieve further 
understanding it is useful to assign the remaining ■0, x fields on the odd sublattice 
A Q certain graphical symbols as follows (The full black dot denotes a point x on the 
even sublattice while the hollow circle stands for the point on the odd sublattice 
which the argument of the fields relates to.). 

o— = i)x,i{x-e x ) (31) o-~ = ^ 2 ,iO-ei) (32) 
•hmd = ^i{x + ei) (33) »-hm3 = 0i )2 (x + ei) (34) 
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I = Xi,i(z-e 2 ) (35) | = X2,i{x-e 2 ) (36) 

J = X2,2(x + e 2 ) (37) J = Xi,2^ + e 2 ) (38) 



As a rule, an arrow flowing out of a point of the odd sublattice (hollow circles) 
denotes a field with flavour index 1 while an arrow flowing into a point on the 
odd sublattice relates to a field with flavour index 2. From this it is already clear 
that only those products of (four) Grassmann fields at any point of the odd sub- 
lattice allow non-vanishing contributions to the partition function Z\ which have 
two incoming and two outgoing arrows in their graphical symbols. We will refer 
to this fact as the vertex arrow rule. The symbols for the fields on the odd sublat- 
tice can now be used in a natural way to construct further graphical building blocks. 



o 



6 



o 



o 



+ 



o 



6 



Xi,i( x - e s) X2,2(x + e 2 ) + XiA x + e 2) X2,i{x - e 2 ) 



(39) 



o — • — o = 



o » • - o 



+ 



O ■ • ' o 



^i,i(a;-ei) ^2,2(3? + ei) + ^1,2(2; + e x ) if} 2 ,i(x - e ± ) 



(40) 





+ 




Xi,i(x - e 2 ) ^2,2(2; + ei) + ^1,2(2: + ei) X2,i(x - e 2 ) 



(41) 





+ 




ipi,i(x - ei) X2,2( x + e 2) + Xi^{x + e 2 ) ip 2 ,i(x - e x ) 



(42) 
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+ 

XiA x + e z) ^2,2{x + ei) + ipi,2(x + ei) X2,2( X + e ?) (43) 

+ 

= ^i,\{x - e x ) X2,i{x - e 2 ) + - e 2 ) ^2,1 - e i) (44) 

1 = I x 1 Xx,i{x - e 2) X2,i(x ~ e 2 ) (45) 

J = I x J = X iA x + e a) X2,2(x + e 2 ) (46) 

o— • = o— • x o-k = tpi,i(x - ei) ^2,1 (a - ei) (47) 

•— = -—-CD x «-^-o = ^i, 2 (x + ei) ^ 2)2 (x + ei) (48) 

Using above pictorial language for the remaining fields and their combinations on 
the odd sublattice we can gain now already some intuitive understanding for the 
structures appearing on the r.h.s. of eqs. (p4])-pB]). We give some characteristic ex- 
amples, similar relations can also be obtained for the vertices omitted now. Please 
note, that the symbols on the l.h.s. below stand for the original expressions (|TTD, 
(HID, (HI), (HID, ( HP restricted to the ip, x factors (The numerical factors have been 
included with their roots, therefore.). 



Vertex 3 



o > « > o x o ■ • ' o = - o— • x —o (49) 



Vertex 6 



x « (50) 
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Vertex 10 



1 
2 




x 




o — • O X 



(51) 



Vertex 14 



V2 



o — • O X 




V2 




(52) 



Vertex 20 



1 

2 




x 






(53) 



One immediately recognizes the significance of the numerical factors for each ver- 
tex. The rule simply is that each corner element (f4"T|)-([44|) is associated with a factor 
l/y/2. Although more useful at the end, in this respect the graphical representations 
given rightmost in eqs. (150), (]5T|) are somewhat misleading. 



What remains to be done now is to discuss the result for the Grassmann in- 
tegration on the odd sublattice. For the moment we exclude vertex 20 (and its 
equivalent on the odd sublattice) from the consideration as its complexity requires 
special attention. It will be discussed at a later stage. In view of the vertex arrow 
rule one finds the following non-zero results (We give only few typical cases, all 
others are related by symmetry considerations of the pictures using rotations and 
reflections.). The graphical symbols below denote in an intuitive way the product 
of (four) Grassmann variables, their ordering performed clockwise with respect to 
the points on the even sublattice starting with the leftmost point in a picture. The 
black dot in the hollow circle of the point on the odd sublattice indicates that the 
Grassmann integration has been performed now. 



2 2 



= / II II dt PiA x ) ^1,2(2) ifa,2(x) ^i,i{x) ^2,1 (a?) = -1 (54) 
L = -l (55) 
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» « I < • — — • » 1 ; » • 



V2 



(56) 
(57) 



From these results one recognizes that links with ki — 1 (thin lines) must form 
self-avoiding loops. It remains to determine the weight of such loops inasmuch as 
vertices contributing to such loops have two terms each (see, e.g., eqs. flSll) , (|52|) 
and (|40|), (^2|)). However, due to the vertex arrow rule only those products of terms 
give a non-vanishing contribution which form an oriented loop. This leaves for each 
loop just two terms (corresponding to the two possible orientations of the loop). 
Both of these terms give exactly the same contribution. This can be seen easily by 
interchanging flavours 1 and 2 what amounts to arrow reversal. Each of the building 
blocks of a thin line (hi = 1) loop (E0)-(H3) transforms under flavour exchange into 



the negative of itself. Also, in view of eqs. flog ), (|57|) at each point belonging to 
the odd sublattice of a given thin line loop arrow reversal brings a relative minus 
sign about. Now, these two minus signs match locally at each point of the odd 
sublattice of a given thin line loop. Consequently, the weight of each thin line loop 
is given by the product of the weights of its constituting vertices multiplied by a 
factor of 2 (the loop multiplicity). Finally, one recognizes from eqs. (|5"5])-(|5lj) that 
each corner (being the product of 2 fields originally) contributes a factor of l/v^2 
(For eq. QoTl| ) one has to apply an interpretation analogous to the middle part of eq. 
(|5TD, of course.). 



Let us now extend the analysis by also allowing configurations containing vertex 
20 we omitted so far. It is clear from the beginning that the picture of self-avoiding 
thin line (hi = 1) loops cannot be maintained in any naive way, but we will see 
how it can still be used in a modified version. The above expression (|53|) already 
suggests what can be done. One may identify the vertex 20 with two loop segments 
passing each other in a given point without intersection. Clearly, there are two ways 
of doing so which are reflected in the two different terms shown in eq. (^). Each of 
these two terms corresponds to four possible combinations of segments of oriented 
loops (each being related to the product of four fields on the odd sublattice, cf. eqs. 
(f4TD-(pBD). If the two segments happen to be just different parts of one loop two of 
the four possible combinations of oriented thin line loop segments are effective in 
the final Grassmann integration over the odd sublattice only (The other two com- 
binations must then lead to a zero result according to the vertex arrow rule.). If 
they belong to two different loops all terms contribute at the end. In considering 
the two terms in eq. fl53|), it is clear however that, if one term is part of two different 
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loops, the other must just belong to one loop returning at some stage to the point 
x G A e under consideration. Consequently, the two terms in eq. ([53]) differ in their 
contribution at the end just by a factor of 2 if it comes to the loop count (just for 
geometrical reasons). As each self-avoiding loop contributes a factor of 2 (related 
to the two different orientations of the loop) the term which is part of two different 
loops will be related to a factor of 4 while the other one is related to a factor of 2 
only. The remaining worry is what can be said with respect to the relative sign of 
the contributions from the two terms in eq. (|53|). We will see that the relative minus 
sign in eq. ( |53"|) is spurious (arising from ordering effects of Grassmann variables) 
and that in fact the contributions from the two terms in eq. add up at the end. 
This can be detected by noting that if one writes the two geometrically different 
terms in eq. ( [53"D as sum over products of (four) fields two of these products are 
the same, namely 4>i t i(x — e i) X2,2( x + ^2) i>i,2{ x + e%) X2,i( x — ^2) (The horizontal 
arrows are flowing in into the point x while the vertical ones are flowing out.) and 
^2,1 s — e i) Xi,z{ x + 62) i/>2,2( x + e i) Xi,i( x ~~ e 2) (The horizontal arrows are flowing 
out of the point x while the vertical ones are flowing in.). Both contribute to the 
l.h.s. of eq. (|53| ) with a total factor of 1 at the end. This proves that contributions 
from both terms in eq. (^) add up after properly taking into account ordering effects. 

Finally, to complete the picture we have to discuss the analogue of vertex 20 
on the odd sublattice. One finds the following non-vanishing Grassmann integrals 
beyond those given in eqs. (|54"|)-(|5"7|) (We apply here the same conventions as there.). 




Again, the loop picture is helpful to recognize the underlying structure. Each of the 
vertices given in eq. (|5^) allows in an unique way an interpretation as standing for 
two segments of oriented loops passing each other in the given point on the odd sub- 
lattice without intersection (The first two diagrams are standing in correspondence 
to the first term in eq. (^) while the other two are being analogous to the second 
term in eq. ([53]).). This interpretation is supported by the appearance of the factor 
1/2 (characteristic for two corners) on the r.h.s. of eq. fl5"9|) . Arrow reversal for just 
one loop corner transforms any diagram of eq. ( |59"D into one of the diagrams shown in 
eq. (|58|). This immediately sheds light onto the different numerical result for these. 
We may imagine each diagram in eq. ( |58P being a sum of two identical terms each 
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of which is supplied with an additional specification. The additional specification is 
provided by one of the two possible ways of two oriented loop segments passing each 
other without intersection in the given lattice site (The point is that the diagrams in 
eq. (|58|) in opposition to those in eq. ( [59"D allow two different interpretations each.). 
Each of those terms with the additional specification contributes —1/2 then. In 
analogy with eq. fl57|), under arrow reversal for just one oriented loop corner we have 
obtained this way a partner for each diagram in eq. (^). This partner diagram has 
as weight the negative of the corresponding diagram in eq. (|59"D. This completes the 
discussion. 

The above consideration allows us to use the following efficient graphical rule to 
determine the weight of any cluster containing vertex 20 (and its analogue on the 
odd sublattice). 




In words, the rule shown above simply requires to replace any vertex 20 (and its 
analogue on the odd sublattice) within a given cluster by two vertices built of two 
thin line loop segments passing each other without intersection. This leads to 2 n 
different cluster configurations when n is the number of times vertex 20 (and its 
analogue on the odd sublattice) occurs in the original configuration. The investiga- 
tion of the ip, x subsystem shows that this rule can reproduce the contribution of 
vertex 20 (and its analogue on the odd sublattice) if one assigns each of the two new 
vertices a factor 1/2, counts any thin line (ki = 1) loop in the newly created con- 
figurations with a factor (loop multiplicity) of 2 (characteristic for oriented loops), 
and then adds up the weights of all these configurations. As the loop count gives 
different results in a given geometrical situation for the two vertices appearing on 
the r.h.s. of eq. (^) the generalization of the above rule to the whole ■?/>, x, ip, \ 
system requires further thought. As the rule fl60|) can be applied for the ■?/>, % an d 
the ip, x subsystems separately each occurrence of a vertex 20 (and its analogue on 
the odd sublattice) leads to a situation where it is met twice in a product of field 
combinations related to this lattice site as one recognizes from the picture below 
(The vertices of the ■0, x subsystem (upper line) are emphasized by stronger lines. 
If one wants to have in mind some explicit expression for the intuitive picture below, 
eq. fl5"5| ) is recommended as reference.). 
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X 




(61) 



There does not exist any problem connected with the two terms where the vertices 
are identical in the -0, x and the ip, x subsystems. However, there are also two cross 
terms which have to be taken care of. In principle, it would be our aim to apply the 
rule (^) to the complete ■0, x, ■0, x system by simply changing the thin line loop 
multiplicity from 2 (for each subsystem) to 4 (for the complete system) but these 
cross terms prevent us from doing so. Consequently, if one wants to determine the 
weight of a given cluster for the full 0, x, -0, x system one has to carry out the loop 
count for one subsystem with loop multiplicity 2 and then to square the result to get 
the correct final weight. The study of some simple clusters suggests that one might 
be able to simplify this procedure by redefining the weights of the two vertices on the 
r.h.s. of eq. ( |60D (i.e., by hiding the contribution of the two cross terms in eq. (|6l|) 
in a modified weight for the two terms where the vertices are identical in the -0, x 
and t/>, x subsystems). However, closer inspection reveals that such a simplification 
is not possible. 

We are now in a position to characterize the vertex model equivalent to the 
two-flavour strong coupling lattice Schwinger model with Wilson fermions. It is a 
(modified) 3-state 20-vertex model on the square lattice. The vertices are shown 
in Fig. 1 and the vertex weights are given below. (All other vertices of the 3-state 
model one might think of have vanishing weight.). The term 3-state refers to the 
three possible states of links (ki = 0, 1, 2; dashed, thin and thick lines respectively). 
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Cl>4 — CU5 — 



u 6 = u 7 = u 8 = u 9 = u 10 = u u = u 2 o = - (64) 
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(65) 
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Figure 1: Vertices of the (modified) 3-state 20- vertex model on the square lattice 
equivalent to the two-flavour strong coupling Schwinger model 



The partition function reads 



Z A [z] 



E 



|A|-|L| 



N(Lf 



l\C(L) 
2J 



(66) 



with 



N(L) = £ N(s,L)2 s . (67) 

s=0 

The sum in eq. (^) is extended over all possible configurations L that can be built 
from the vertices 1-20. |A| is the number of lattice points and \L\ the number of 
links with thin or thick lines (ki = 1, 2). C(L) counts the number of times a thin line 
bends (For the purpose of this count thick lines are understood as being equivalent 
to two thin lines and the vertices 8-11 therefore represent a twofold bending of a 
thin line. For vertex 20, we recall the rule flB0|). P n {L) is the number of times the 
vertex n appears in a given configuration L.). 

11 19 

C{L) = 2 £ P n {L) + ]T P n {L) + 2 P 20 (L) (68) 

n=4 n=12 
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The factor 1/2 then can be understood as a bending rigidity of thin lines. The same 
factor for the bending rigidity has been found in the one-flavour strong coupling 
Schwinger model N(s, L) is the number of configurations with s thin line loops 
generated from L by application of the rule fl60|). Following sum rule of course holds. 

£ N(s,L) = 2 P2o(i) (69) 

s=0 

If P2o(L) = 0, one finds therefore 



N{Lf = (70) 

where S(L) denotes the number of thin line loops in a given configuration L. In 
this simple case, thin line loops enter with a multiplicity of 4 (This case bears some 
similarity to the general loop model discussed in 0.). 



Eq. fl66|) defines a modified 3-state model because in the usual definition of a 
3-state model the factor N(L) 2 would not be present. Q-state models with certain 
similarities to eq. ( |66"D (without the factor N(L) 2 ) have been investigated e.g. in 
[P]-[|18|. For z = we obtain a 19- vertex model. Somewhat similar 19- vertex models 
(without the factor N(L) 2 ) have recently been dicussed in [19|-[2~I]. Now, the mod- 



ified 3-state model flBTp can alternatively be understood as a regular 6-state model. 
To see this note, that with reference to the graphical rules given above each thin line 



(ki = 1) can be considered as having 4 different states. According to eqs. (pl|)- (|38|) 
each thin line term considered within the ■?/>, \ and ip, \ subsystems can be related to 
an arrow assigned to this line. This assignment leads to four different states for thin 
lines and transforms the 3-state model into a 6-state model. Consequently, in this 
picture the vertices 8-19 split into four different 6-state model vertices (each corre- 
sponds to one of the four different products of fields related to a given 3-state model 
vertex, cf. eqs. (]T6|) (p7[) ) • Each of these vertices has the weight of the corresponding 
3-state model vertex. The vertex 20 yields 6 x 6 = 36 different 6-state model vertices 
(in each subsystem the vertex 20 corresponds to 6 different products of fields only, 
cf. eq. (|28|) ). The weights of the 6-state model vertices derived from the vertex 20 
can easily be determined by expanding the r.h.s. of eq. ( p8|) in terms of products 
of fields. Taking into account the above considerations one may conclude that the 
modified 3-state 20-vertex model fl6"6|) is equivalent to a regular 6-state 91-vertex 
model. To which extent this is actually useful information remains to be seen in the 
future. Another line of thought might be linked to the following consideration. As 
N(L) is related to some loop count the vertex model (|66D bears some resemblance to 
a loop model. The loop model picture becomes even more obvious if one interprets 
thick lines in the present model as two parallel thin lines. If one then applies the 
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rule ( pOD to both the ijj, \ and the ip, x subsystems and assigns each subsystem a 
colour, we end up with a 2-colour loop system on the square lattice. A loop model 
which comes closest to this picture has been discussed in P^] . However, the relation 
to the present model does not go beyond a certain remote similarity. 

Finally, in analogy with the one-flavour case let us discuss the critical behaviour 
of the two-flavour strong coupling Schwinger model by means of some approximate 
method. We are going to rely on a generalization of the independent loop ap- 
proximation (cf. sect. 3 of 0) which we now appropriately call independent cluster 
approximation. As in the case of the independent loop approximation we assume 
that the system under consideration is sufficiently dilute in the relevant parameter 
region in order to yield reasonable results. We start with the observation that any 
cluster built of the vertices 2-20 falls into two classes. Either we have a self-avoiding 
thick line loop or some cluster containing thin line segments. In the independent 
cluster approximation we now write the partition function Z&[z] as follows. 



ZaW = z ,m e {^[W/4] + Z A [1,« 



(71) 



Here, Z&[1, z, r/] is the single loop partition function defined in sect. 3 of 0. Z\[l, z] 
denotes the single cluster partition function in which the statistical sum is extended 
only over all possible configurations L containing just one cluster containing thin 
lines. z' A ' Z\[l, z, 1/4] coincides with Z\[z] when the summation in eq. ( |66|) is 
restricted to configurations L with just one thick line loop. z' A ' Za[1,z) coincides 
with Zjs\z\ when the summation in eq. (|66|) is restricted to configurations L which 
constitute just one cluster with thin line elements. As the free energy density / 
reads in the independent cluster approximation {(3t = 1/T) 

Prf(z) = - lim -t- \nZ A [z] 

|A| >oo |A| 

= _,„,_ , im g^VjL+AM (72) 

|A|— ^oo |A| 

a phase transition can be detected by looking at each single cluster (loop) partition 
function separately. We do not have any information about Z\[l,z] but the single 
loop partition function Zx[l,z,r]] has been studied in ||. One finds a critical line 
described by the equation 



z(r)) = 1 + ^-1) . (73) 
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/i = 2.638... is the effective coordination number of the square lattice. From eq. ([73]) 
one immediately obtains in the present case 



M\M\ 



1 



Gu + 3) 



(74) 



■cr 



I6K1K2 
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From the experience with the one-flavour JT| and the two-flavour strong coupling 
Schwinger models it seems reasonable to expect that for the investigation of the 
critical behaviour of the general Ay-flavour model the self-avoiding loop model with 
bending rigidity rj = 2~ N f will be relevant. 

To conclude, we have found the vertex model equivalent to the two-flavour strong 
(infinite) coupling lattice Schwinger model with Wilson fermions. This demonstrates 
that the method employed by Salmhofer for the one-flavour case |l| can success- 
fully be generalized (with some technical effort) to models with a higher number of 
Grassmann variables per lattice site (8 in the present case; incidentally note that 
one-flavour strong coupling lattice QED with Wilson fermions in 4D has the same 
number of Grassmann variables per lattice site). As the method can also be applied 
to purely fermionic models this suggests a number of interesting directions for fur- 
ther research which will be explored in the future. 

Acknowledgements 

The present work has been performed under the EC Human Capital and Mo- 
bility Program, contract no. ERBCHBGCT930470. I would like to thank Simon 
Hands for discussions and a critical reading of the draft version of the paper, and 
Holger Perlt for valuable advice concerning Mathematica. 



20 



References 

M. Salmhofer; Nucl. Phys. B 362(1991)641. 
P. Rossi, U. Wolff; Nucl. Phys. B 248(1984)105. 
H. Gausterer, C.B. Lang; Phys. Lett. B 341(1994)46. 
F. Karsch, E. Meggiolaro, L. Turko; Phys. Rev. D 51(1995)6417. 
H. Gausterer, C.B. Lang; Nucl. Phys. B 455[FS] (1995)785. 
K. Scharnhorst; University of Wales, Swansea, Preprint SWAT/95/72, |hep- 



[2] 
[3] 
[4] 
[5] 



[6 



[7; 



[9 



[at 9505001 



Wolfram Research; Mathematica, version 2.2. Wolfram Research, Champaign, 
1994. 

F.S. Rys, W. Helfrich; J. Phys. A 15(1982)599. 
Yu.G. Stroganov; Phys. Lett. A 74(1979)116. 



[10] J.H.H. Perk, C.L. Schultz; Phys. Lett. A 84(1981)407. 

[11] J.H.H. Perk, C.L. Schultz, in: M. Jimbo, T. Miwa (Eds.): Non-linear Integrable 
Systems — Classical Theory and Quantum Theory, Proc. of RIMS Symposium 
organised by M. Sato, Kyoto, Japan, 13-16 May 1981. World Scientific, Singa- 
pore, 1983, p. 135. 

[12] K. Sogo, Y. Akutsu, T. Abe; Prog. Theor. Phys. 70(1983)730, ibid. 739. 

[13] J.H.H. Perk, C.L. Schultz; Physica A 122(1983)50. 

[14] J.H.H. Perk, F.Y. Wu; Physica A 138(1986)100. 

[15] H.J. Giacomini; Phys. Lett. A 117(1986)381. 

[16] Y. Akutsu, M. Wadati; J. Phys. Soc. Jpn. 56(1987)3039. 

[17] L.H. Gwa, F.Y. Wu; J. Phys. A 24(1991)L503. 

[18] M.T. Batchelor, J. Suzuki; J. Phys. A 26(1993)L729. 

[19] M. Idzumi, T. Tokihiro, M. Arai; J. Phys. I France 4(1994)1151. 

[20] Y.M.M. Knops, B. Nienhuis, H.J.F. Knops, H.W.J. Blote; Phys. Rev. B 
50(1994)1061. 



21 



[21] A. Kliimper, S.I. Matveenko, J. Zittartz; Z. Phys. B 96(1995)401. 
[22] S.O. Warnaar, B. Nienhuis; J. Phys. A 26(1993)2301. 



22 



